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Abstract. We determine the Strouhal number (hereafter St), which is essentially a nondimensional measure of the correlation 
time, from numerical calculations of convection. We use two independent methods to estimate St. Firstly, we apply the 
minimal tau-approximation (MTA) on the equation of the time derivative of the Reynolds stress. A relaxation time is obtained 
from which St can be estimated by normalising with a typical turnover time. Secondly, we calculate the correlation and 
turnover times separately, the former from the autocorrelation of velocity and the latter by following test particles embedded 
in the flow. We find that the Strouhal number is in general of the order of 0.1 to 1, i.e. rather large in comparison to the typical 
assumption in the mean-field theories that St -C 1. However, there is a clear decreasing trend as function of the Rayleigh 
number and increasing rotation. Furthermore, for the present range of parameters the decrease of St does not show signs of 
saturation, indicating that in stellar convection zones, where the Rayleigh numbers are much larger, the Strouhal number may 
indeed be significantly smaller. 
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1. Introduction 

The mean-field theories of hydromagnetic dynamos and 
angular momentum transport need knowledge of turbulent 
correlations, namely the electromotive force and Reynolds 
stresses, respectively. The explicit calculation of these is vir- 
tually impossible for any astrophysical conditions due to the 
lack of necessary computational cababilities. Therefore the 
small-scale effects are usually parametrised by transport co- 
efficients which relate the turbulent correlations to the mean 
quantities (e.g. Steenbeck, Krause & Radler 1969 Riidiger 
19891. However, in order to calculate the transport coeffi- 
cients, knowledge of the small-scale velocities and magnetic 
fields are still needed. Thus, simplifying assumptions such 
as the first order smoothing approximation (hereafter FOSA) 
have been used (Steenbeck et al. 1969 1. Although the results 
obtained with FOSA are in agreement with observations in 
many cases, the basic assumptions of the theory have rarely 
been thoroughly studied (however, see Petrovay & Zsargo 
1998 Kapvla et al. l2004b> . 
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The validity condition for FOSA is that either the 
Reynolds or the Strouhal number is small. The former con- 
dition is clearly never met, but of the latter in the context of 
convection and dynamo theory no systematic study seems to 
exist. A further notion associated with the Strouhal number is 
that even though if St is not small in the sense that St <C 1, 
it is still possible to calculate the transport coefficients from 
a cumulative series expansion if St < 1 (e.g. Knobloch l 19781 
Nicklaus & Stix fT988b . 

Recently it has been shown numerically that St can ex- 
ceed unity for forced turbulence (Brandenburg, Kapyla & 
Mohammed 2004 Brandenburg & Subramanian 2004) in the 
contexts of passive scalar diffusion and mean-field dynamos, 
respectively. In these studies, instead of FOSA, the minimal 
tau-approximation (hereafter MTA) is found to be in better 
agreement with the calculations. In MTA, instead of the tur- 
bulent correlation itself, the time derivative is investigated, 
and the higher order terms are parameterised by a term which 
is just the original turbulent correlation divided by a relax- 
ation time (Blackman & Field 2003 1. Furthermore, interpret- 



2 THE STROUHAL NUMBER 



2.2 Determination of the timescales 



ing the relaxation time as the correlation time of the turbu- 
lence the Strouhal number can be estimated. 

The forced turbulence results raise the question of the 
value of St for convection, which seems to be rather badly 
known. For example, solar surface observations indicate that 
the lifetime and turnover time of granules is approximately 
the same, yielding St w 1 (e.g. Stix 2002 1. However, this re- 
sult may not be relevant for the solar dynamo which is work- 
ing in the deeper layers. In the present study we estimate the 
Strouhal number for numerical convection by two indepen- 
dent methods. Firstly we apply the MTA on the equation of 
the Reynolds stresses, and secondly we calculate the corre- 
lation and turnover times separately directly from the flow. 
The correlation time is estimated from the autocorrelation of 
velocity and the turnover time by following test particles em- 
bedded in the flow. The computational model is a Cartesian 
box situated at a latitude on a star. The model is described 
in detail in Kapyla, Korpi & Tuominen ( 2004a I. 

The remainder of the paper is organised as follows: in 
Sect. (0 the two methods used to estimate the Strouhal num- 
ber are discussed. Sects. (0 and (|4j give the results and con- 
clusions, respectively. 

2. The Strouhal number 
2.1. Minimal tau-approximation 

We apply the MTA on the time derivative of the Reynolds 
stress, Qij = (u^u'j), where the brackets denote horizontal 
averaging and primes the fluctuation. We arrive at the equa- 
tion (see Kapyla et al. 2004b for a more detailed derivation) 



dQij 
dt 
where 



Trel 



(1) 



9 ijk = -2 (cikiiu'^) + e iW («{«{» , (2) 

where ejjy is the Levi-Civita symbol and the rotation vec- 
tor. Taking into account only the lowest order effect for the 
vertical A-effect, which is proportional to the component Q yz 
in the local convection model, the relaxation time turns out to 
be 



2£lcosO(Q zz - Qyy) 

Interpreting T rc i as the correlation time of the turbulence, the 
Strouhal number can be calculated identically as in the forced 



turbulence studies (Brandenburg et al. 2004, Brandenburg & 
Subramanian 2004 1 

St^ = kfU Tms T re i , (4) 

where kt is the wavenumber of the energy carrying scale and 
M rms the average rms-velocity in the convectively unstable re- 
gion. For convection calculations k{ corresponds essentially 
to the largest possible scale permitted by the box dimensions 
(see Kapyla et al. l2004bt 
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Fig. 1. The Strouhal number from MTA. 



2.2. Determination of the timescales 

The disadvantage of the MTA-approach is that it can only be 
applied to the case where rotation and the Reynolds stresses 
are statistically nonzero. In order to circumvent this problem, 
we have devised an independent way to estimate St which 
works also if rotation is not present and without any recourse 
to the mean-field theory. The method essentially consists of 
the calculation of the correlation and turnover times sepa- 
rately and directly from the flow. 

The correlation time is estimated from the autocorrelation 
of velocity 



C[ui(x,t ),Ui(x,t)} 



Ui(x,t )Uj(x,t) 
y/u?(x,t Q )uj(x,t) 



(5) 



where x is the position vector and to and t denote the times 
from which the snapshots were taken. The correlation time r c 
is defined as the time after which the correlation drops below 
a threshold value. In the present study the threshold is set to 
0.5. We use the vertical velocity component to determine the 
correlation time. 

The turnover time, t% Q , is estimated from the trajectories 
of test particles which are advected by the flow. The turnover 
time can be defined as the time which passes between two 
consecutive crossovers of some fixed reference level into the 
same direction. However, this definition implicitly assumes 
that the vertical scale of convective motions is of the order of 
the depth of the convectively unstable layer. Thus, turnovers 
happening far away in comparison to the scale of convection 
are not registered at all. Another possibility is to define the 
turnover time as the time which elapses between two con- 
secutive changes of direction (into the same direction). This 
latter definition registeres all turnovers and we shall present 
results using that in the remainder of the paper. Differences 
between the results obtained with the two definitions are dis- 
cussed further in Kapyla et al. (2004b). 

Once the correlation and turnover times have been deter- 
mmined, the Strouhal number is simply their ratio 
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Fig. 2. The Strouhal number as function of the Rayleigh (left) and Coriolis (right) numbers. The dashed lines give the power 
laws St oc Ra" 49 and St oc Co" 45 



3. Results 

The minimal tau-approximation can be applied to the 
Reynolds stresses only if there is appreciable rotation, and 
away from the poles where the stresses vanish due to sym- 
metry in the present geometry (e.g. Kapyla et al. 2004a I. A 
typical result is shown in Figure[T] The analysed calculations 
are those made at a latitude 30 degrees south from Kapyla et 
al. J200 4ai. The stress is taken from the middle of the convec- 
tively unstable layer which generally quite well describes the 
situation in the whole box. The errorbars denote the modified 
mean error of the stress (see Eq. (30) of Kapyla et al. 2004a). 
As function of the Coriolis number the Strouhal number is 
essentially determined by the value of r re i, which varies in a 
very similar manner as the vertical A-effect. The reason for 
this behaviour is that u rms and the diagonal Reynolds stresses 
Qyy and Qzz vary only little as function of Co leading to the 
fact that the functional form for the relaxation time and the 
A-effect is essentially the same, Q yz /Q. The values of St( rcl ) 
vary from about 0.4 for Co = 2 to essentially zero for very 
slow and very rapid rotation. 

Whereas the applicability of the MTA-approach is limited 
to cases with rotation, the Strouhal number should have a fi- 
nite value also without rotation. Thus, we set out to extract 
the correlation and turnover times separately from the flow as 
described in Sect. i2.2\ . We find that the the correlation time 
decreases consistent with r c oc Ra -0 ' 45 as function of the 
Rayleigh number. The turnover time, however, changes only 
marginally. Thus the Strouhal number follows approximately 
a power law St oc Ra (see left panel of Figure [2}. Al- 
though the parameter range in the present study is quite lim- 
ited, the result is still promising in the sense that if the same 
trend carries over to stellar parameters, the Strouhal number 
may be much smaller there. 

A similar trend is found when rotation is increased. The 
correlation time decreases rapidly as function of the Cori- 
olis number whereas the changes in the turnover time are 
only minor. However, one aspect not to be overlooked here is 
the trend seen in the turnover time. The simple estimate, the 



depth of the convectively unstable region divided by the aver- 
age velocity, increases as function of rotation due to the fact 
that overall velocities tend to diminish as rotation becomes 
more rapid. However, the turnover time calculated from the 
test particle trajectories shows an opposite trend due to the 
fact that the spatial scale of convection is reduced even more 
than the overall velocities. Thus we find that for moderate and 
rapid rotation, the Strouhal number approximately follows a 
power law St oc Ra -0 ' 5 (see right panel of Figure [2). The 
explanation is most probably the strong Coriolis forces which 
tend to disrupt any coherent flow structures. 

4. Conclusions 

We estimate the Strouhal number from numerical models of 
convection with two independent methods. Firstly, we ap- 
ply the minimal tau-approximation of the equation of the 
Reynolds stress. Secondly, we calculate the correlation and 
turnover times directly from the flow without any recourse to 
a mean-field theory. 

We find that the Strouhal number from the MTA reaches 
values of maximally w 0.5. St has a maximum for intermedi- 
ate rotation, where the Reynolds stress itself also peaks. The 
Strouhal number follows closely the same trend as the ver- 
tical A-effect (see e.g. Kapyla et al. 2004a) as function of 
rotation due to the similar functional form. 

The correlation time is seen to decrease consistent with 
power law t c oc Ra~°' 4j as function of the Rayleigh num- 
ber and approximately with r c oc Co~° 45 for moderate and 
rapid rotation. These results indicate that although the val- 
ues of St in this study are generally of the order of 0.1 to 1, 
the Strouhal number in stellar convection zones, where Ra is 
most certainly, and Co at least probably, much larger than in 
the present study, may be significantly smaller. 
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